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Influence of an inter-chain coupling on spiral ground-state correlations in frustrated 
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We investigate the influence of an inter-chain coupling on the spiral ground state correlations of 
a spin- 1/2 Heisenberg model consisting of a two-dimensional array of coupled chains with nearest 
(Ji) and frustrating next-nearest neighbor (J2) in-chain exchange couplings. Using the coupled- 
cluster method we calculate the transition point between the commensurate and the incommensurate 
(spiral) ground state as well as the quantum pitch angle of the spiral ground state. In addition, we 
provide a simple empirical formula which describes the relation between the quantum pitch angle 
and the frustration parameter J2/J1. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.40.Mg, 75.50.Eo 



I. INTRODUCTION 

In recent years strongly frustrated quantum magnets 
exhibiting exotic ground state (GS) phases have been 
intensively investigated both theoretically and experi- 
mentally, see e.g. Refs. [l|, [3, S At zero temperature 
all transitions between GS phases are driven purely by 
the interplay of quantum fluctuations and the compe- 
tition between interactions (e.g., frustration), see e.g. 
Refs. 3, y, [y, 0, & Particular attention has been paid 
to one-dimensional (ID) J1-J2 quantum Heisenberg mod- 
els, which have been studied theoreticallywith much suc- 
cess over the last two decades, see Ref. and references 
therein. From the experimental side recent investigations 
have shown that edge-shared chain cuprates build a spe- 
cial family of frustrated quantum magnets which can be 
described by a quasi-lD J1-J2 Heisenberg model. Among 
others we mention here LiVCu04, LiCu202i NaCu202, 
Li2ZrCu04, andLiaCuOs 0, 13, [III, [H, El, HE HE 
d,llMl> which were identified as quasi-lD frustrated spin- 
1/2 magnets with ferromagnetic nearest-neighbor (NN) 
in-chain Ji and antiferromagnetic next-nearest-neighbor 
(NNN) in-chain interactions J2. These cuprates have at- 
tracted much attention due to strong quantum effects 
and the observation of incommensurate spiral (helical) 
spin-spin correlations at low temperatures. Among these 
materials Li2ZrCu04 and Li2Cu02 are of particular in- 
terest, since these compounds are found to be near a 
quantum critical point [iTI . Il9l | . 

The ID frustrated spin-1/2 J1-J2 Heisenberg model 
may serve as the simplest model to describe some im- 
portant features of such materials. The GS properties of 
the model in the classical limit are well known, i.e. when 
the spin quantum number s — » 00. In this case the GS 
does exhibit a second-order transition from a coUinear 
phase (ferro or antiferro) to a non-coUinear phase with 
spiral correlations along the chains at J2 = |Ji|/4. For 
•h > 1^1 1/4 the classical spiral (pitch) angle a^ is given 

by 



ckci = arccos( — J1/4J2) . 



(1) 



Similar expressions can be derived in the presence of fur- 
ther couplings J3, J4, etc., however, this more general 
case will not be considered here for the sake of simplic- 
ity. Note that in the classical limit neither the pitch 
angle ad nor the transition point J2 = | Ji|/4 depend on 
the inter-chain coupling J^. Such a classical relation as 
Eq. ([T]) between Ji , J2 , ... and the pitch angle aci has 
been used to justify ID parameter sets obtained from 
fitting x{T) data [20| or by LDA mapping procedures 
[2l| . However, ignoring this way sizable quantum effects, 
such an approach is not very convincing. Furthermore, in 
real materials additional terms in the Hamiltonian such 
as anisotropy or exchange coupling between the chains 
might be of relevance to yield a quantitative theoretical 
description of the experimental results. The existence of 
helical long-range order at low temperatures makes the 
importance of the inter-chain coupling evident. However, 
from the theoretical side such extended models so far are 
much less studied than the pure ID-" parent" models. 

Therefore, in the present paper we focus on the discus- 
sion of the effect of the inter-chain coupling J±^ on the 
GS spin-spin correlations in frustrated spin-1/2 Heisen- 
berg model consisting of a two-dimensional (2D) array 
of frustrated chains coupled by J_l. In particular, we 
will discuss GS's with incommensurate spiral (i.e. non- 
collinear) correlations. Though, meanwhile many papers 
exist dealing with spiral correlations in the strictly ID 
system (J_l = 0) [23,[2l,[2l,[25, 26, 27, 28], the influence 
of the inter-chain coupling on the pitch angle to the best 
of our knowledge has not been discussed so far. 

The theoretical treatment of frustrated quantum an- 
tiferromagnets is far from being trivial. Although one 
can find exact GS's in some exceptional cases, see e.g. 
Refs. [22, [2^, [20, [31I, [33, standard many-body meth- 
ods may fail or become computational infeasible. For 
instance, the quantum Monte Carlo techniques suffer 
from the minus-sign problem in frustrated systems. The 
density-matrix renormalization group (DMRG) success- 
fully used to discuss spiral correlations in ID magnets 
[23I . V2d{ is essentially restricted to ID systems, at least 
in the present state of the art. Also the exact diagonal- 



ization technique used in Ref. |26| to find the pitch angle 
for the ID problem would be limited to extremely small 
chain lengths when a finite inter-chain coupling should 
be considered. A method fulfilling the requirement to be 
able to deal with frustrated spin systems at any dimen- 
sion, including magnetic systems with incommensurate 
spiral GS's, is the coupled cluster method (CCM). This 
method was already used for the strictly ID J1-J2 Heisen- 
berg model and it was shown that the CCM results are 
in good agreement with the DMRG data [2^. Hence, in 
the present paper we use the CCM following the lines of 
Ref. [23I but extend the CCM calculations by including 
the inter-chain coupling Jj^. 

The frustrated spatially anisotropic 2D Ji- J2- J_i spin- 
1/2-Heisenberg model considered here reads 



n i 



(2) 



where the index n labels the chains and i the lattice sites 
within a chain n. The NN in-chain coupling Ji is fixed to 
either Ji = 1 (antiferromagnetic) or Ji = — 1 (ferromag- 
netic). The inter-chain coupling J± and the frustrating 
NNN in-chain coupling J2 arc considered as varying pa- 
rameters of the model. Note that in the case considered 
here the inter-chain coupling docs not lead to frustration 
and it is practically of arbitrary strength, i.e. including 
also the region beyond the quasi- ID limit | J^| ^ | Ji|, J2- 
Anyhow, the opposite limit |J^| ^ |=/i|,^2 will not be 
considered for reasons of lacking physical relevance (to 
the best of our knowledge). We focus on the extreme 
quantum case, i.e. (sj „)^ = s{s + 1) with s = 1/2. The 
main point which will be considered here is the influence 
of the inter-chain coupling J± on the transition point 
between the collinear phase and the non-coUinear spiral 
phase and on the pitch angle characterizing the spiral 
correlations in the quantum model with s ~ 1/2. 

II. THE COUPLED CLUSTER METHOD (CCM) 

In this section we outline only some main features of 
the CCM which are relevant for the model under con- 
sideration. Again we mention that we follow the lines 
described in Ref. [23 where the CCM was applied to the 
strictly ID problem. For more details of the method 
the interested reader is referred to Refs. [23, [33, [SJ, [S^, 
m, M, m, M, \M, \M, \M, S, B Special attention to 
the CCM treatment of non-coUinear GS's was paid in 

Refs. M M, M M M \M- 

The starting point for the CCM calculation is the 
choice of a normalized reference or model state 1$), to- 
gether with a complete set of (mutually commuting) 
multi-configurational creation operators {C]^} and the 
corresponding set of their Hermitian adjoints {Cl}, 

($|C+ = = Cl|$) VL^O, Co+ = 1 (3) 

[C+,C+]=0=[Cl,Ck]. (4) 



With the set {|<i>), C^} the CCM paramctrization of the 
ket GS eigenvector |^) of the considered many-body sys- 
tem is then given by 



L#0 



(5) 



The CCM correlation operator S contains the correlation 
coefficients a^ which can be determined by the so-called 
set of ket equations 



mCLe-''He''\<P) =0 ; VL ^ 0. 



(6) 



For the considered frustrated spin system we choose a 
reference state corresponding to the classical state of the 
model, i.e. the ferromagnetic state | mj •••) along 
a chain for Ji = — 1 and small J2 and the Neel state 
I it it • • • ) for Ji = 1 and small J2, whereas for larger 
frustration J2 we have to choose a non-collinear refer- 
ence state with spiral correlations along the chains (i.e., 
pictorially, |$) = | t/'— *\ii/ ■■■)) characterized by 
a pitch angle a, i.e. |$) = |<i>(a)). Such states in- 
clude the ferromagnetic state (a = 0) as well as the 
Neel state {a = n). In the quantum model the pitch 
angle may be different from the corresponding classical 
value Qfci. Hence, we do not choose the classical result 
for the pitch angle but, rather, we consider a as a free 
parameter in the CCM calculation, which has to be de- 
termined by minimization of the GS energy given in the 
CCM formalism by E{a) ^ {<^{a)\e^^He^\^{a)), i.e. 
from dE/da\a=a = 0, see also Appendix, Eqs. (Al)- 
(A6). 

In order to find an appropriate set of creation operators 
it is convenient to perform a rotation of the local axes of 
each of the spins, such that all spins in the reference 
state align in the negative z-direction. This rotation by 
an appropriate local angle Si^n of the spin on lattice site 
(i,n) is equivalent to the spin-operator transformation 



cos Sj 



sin(5i,„sf„; 



(7) 



The local rotation angle (Si.„ can be easily expressed by 
the pitch angle a of the spiral reference state, where the 
relation (5i_„(a) depends on the signs of Ji, J_l and the 
lattice vector Ri^„. 

In this new set of local spin coordinates the reference 
state and the corresponding creation operators C^ are 
given by 



'i,Ti ' °i,n''j,m 1 ^i,n^j,m^k,l ' 



l$) = Uiii 

(8) 

where the indices (i^n),{j,m){k,l), . . . denote arbitrary 
lattice sites. In the new rotated coordinate frame the 
Hamiltonian ^ becomes dependent on the pitch angle 
a. It reads 



Jl 



H = -i ^[cos(a) + l](s+„s,+i_„ + s;;,„s+ 1 „) + [cos(a) - l]is+jt+i,n + h,nh+i,n) + 2sin(a)[s+„sf+i^„ 



s^ s+ 



,n'^i+l,n "^ '^i,n'^i+l,n ^i,n^i+l,n\ + ^ COS(ajSj „Sj_|_]^ , 



■ ^ [cos(2a) + 1] {stnh+2.n + Si.jt+2,n) + [cos(2a) - 1] (s+„s+^2,n + s^.J^+2,n) + 2 sin(2a) [s+„s ■+2,„ 



sf „sio ., + s, „sL, „ - sf ^Sj^2,«] + 4cos(2a)s| „s| 



^i,n''i+2,ri ^ ''i,?i''i+2,n ''i,n'- 



+2,n 



J 



J 



y / X^i.n^i.n+1 + ^i,n^i,n+l + ^Sj „«,; „_[_]^ j + ^^^(Si^nSi^n-j-i + 'Si,n'5i,n+1 + ^^i,n^i,n+l)^ 



(9) 
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FIG. 1: The quantum pitch angle Qqu versus J2 for antifer- 
romagnetic J\ and J± (case (i)). The main panel shows the 
difference between the quantum pitch angle Oqu and its classi- 
cal counterpart ad for various inter-chain couplings J^. The 
inset shows the corresponding data for the quantum pitch an- 
gle Oqu itself. For comparison the classical pitch angle and 



the DMRG data of White and Affleck [25 
quantum model are also shown. 



for the strictly ID 



FIG. 2: The quantum pitch angle Qqu versus J2 for antifer- 
romagnetic J\ and ferromagnetic Jx (case (ii)). The main 
panel shows the difference between the quantum pitch angle 
Oqu and its classical counterpart ad for various inter-chain 
couplings Jj_ • The inset shows the corresponding data for the 
quantum pitch angle aqu itself. For comparison the classical 
pitch angle is also shown. 



where sf,^ = sf_„±isf „, and a {= Si^„~Si+i^„) is the pitch 
angle between the two neighboring spins in a chain inter- 
acting via the NN bond Ji , which has to be determined 
for the quantum model. For ferromagnetically coupled 
spin chains {J^ ^ 0) one has to set J^ = in Eq. ^ 
and vice versa. Therefore, from Eq. ^ it is obvious that 
in the quantum case considered here the role of the (un- 
frustrated) inter-chain coupling J± introduced in Eq. ([2]) 
is different for ferromagnetic and antifcrromagnetic Jj_. 
By contrast, in the classical case a corresponding inter- 
chain coupling does not affect the pitch angle at all. 

The CCM formalism would be exact if we could take 
into account all possible multispin configurations in the 
correlation operator S which is impossible to do in prac- 
tice for a quantum many-body system. Hence, it is nec- 
essary to truncate the expansions of S. In [23| it was 
demonstrated, that the so-called SUB2-3 approximation 



for the strictly ID system leads to results of comparable 
accuracy to those obtained using the DMRG method. In 
the SUB2-3 approximation all configurations are included 
which span a range of no more than 3 contiguous sites 
and contain only 2 or fewer spins. A particular advantage 
of the SUB2-3 approximation consists of the possibility 
to find the relevant CCM equations ^ in closed ana- 
lytical form, see Appendix |^ These explicit equations 
provided here can be used to find the quantum pitch an- 
gle aqu for an arbitrary set of parameters Ji , J2, J± by 
simple numerical solution of them. 



III. RESULTS 

For the model under consideration we have calculated 
that point J2 where the GS state spin-spin correlations 
change from coUinear to non-coUincar spiral correlations 
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FIG. 3: The transition point J2 as a function of the inter- 
chain coupling Jj_ for Ji = 1 and J_l > (case (i)) as well as 
J± < (case (ii)). The triangles indicate a discontinuous and 
the sqaures a continuous change of a^u- 
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FIG. 4: The GS energy versus pitch angle a for different 
(fixed) J2 and antiferromagnetic Ji (Jx =0). 
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FIG. 5: The quantum pitch angle Oqu versus J± for antiferro- 
magnetic Ji = 1 (case (i)) and various values of the frustrat- 
ing NNN exchange J2 (the lines are guides for the eyes) . The 
inset shows the data for J2 — 0.7, 0.8, 0.9 with an enlarged 
scale. 



as well as the quantum pitch angle aqu. In what follows 
we call the point J| the 'transition point'. We mention, 
however, that the question for magnetic GS long-range 
order goes beyond the scope of the present paper. Gener- 
ally one can argue that for the strictly ID problem the GS 
(except the simple ferromagnetic state) does not exhibit 
magnetic long-range order, whereas for finite Jj_ GS long- 
range order can exist, cf. e.g. Refs. |4^, |4^, |43, H, H- 
In particular, the GS phase for larger J2 and small J± 
is magnetically disordered and may have a weak spon- 
taneous dimerization along with finite-range incommen- 
surate magnetic correlations [4^. Hence the 'transition 
point' J2 for small J± may locate that narrow para- 
meter region where finite-range magnetic correlations are 
changing from commensurate to incommensurate ones, 
but does not indicate a true quantum phase transition. 

We present data for J± = ±0, 0.1, 0.2, . . . , 1.0 and a 
fine net of J2 values. For the sake of clarity in the fol- 
lowing we asort the results into four cases, depending on 
the signs of Ji and J±, namely (i) Ji ~ I, J± > 0; (ii) 
Ji = 1, J± < 0; (iii) Ji = -1, J± > 0; (iv) Ji = -1, 
J± < 0. 

Since for cases (i) and (ii) the behavior is quite simi- 
lar, we can discuss both cases together. The pitch angle 
in dependence on J2 is shown in Figs. [T] (case (i)) and 
[2] (case (ii)). For comparison we also draw the corre- 
sponding DMRG data of White and AfHeck [Hi for the 
strictly ID problem in Fig. [TJ These data agree quite 
well with the CCM data in particular for larger J2. It 
is obvious that quantum fluctuations change the classi- 
cal correlations drastically. In particular, by contrast to 
the classical case, the collinear quantum state can sur- 
vive into the region J2 > J2''^ = 0.25, where classically 
it is already unstable. This effect is known as order from 
disorder p'A 5311 and is widely observed in quantum spin 
systems ill Ell [H HI . We find, e.g. J^ « 0.68 for the 
quantum model with J± = 0, which is in good agreement 
with known results, see, e.g.. Fig. 3 in Ref . l26l. 

Switching on J±, i.e. increasing the dimension of the 
spin system, the effect of quantum fluctuations should be- 
come weaker and, as a result, the collinear quantum state 
gives way for the state with incommensurate correlations 
at smaller values of J2, i.e. the transition point moves 
towards the classical value, see Fig. [3l However, even 
in nearly isotropic 2D systems {J± ~ 1) the quantum 
fluctuations are still important and one has J2 ~ 0.38 
(J2 w 0.40) for case (i) (for case (ii)) which is signifi- 
cantly above the classical value. 

The shift of the transition point leads to an interesting 
behaviour of the difference a^^ — ad. It is positive for 
0.25 < J2 < J| but negative for J2 > J^- Surprisingly we 
0.7 as almost independent of J±, i.e. all curves 



find J* 



J^qu 



— ad vs. J2 meet approximately in one point, see 
Figs. [Hand [21 Note that for large J2 > Ji the model with 
strongest quantum fluctuations, i.e. J± — 0, approaches 
the limit aqu -^- 7r/2 most rapidly. 

Next we discuss the behavior of the pitch angle at the 
transition point. For smaller values of the inter-chain 
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FIG. 6; The quantum pitch angle Oqu versus J2 for ferro- 
magnetic Ji and antiferromagnetic J± (case (iii)). The main 
panel shows the difference between the quantum pitch angle 
Qfqu and its classical counterpart a^i for various inter-chain 
couplings J± . The inset shows the corresponding data for the 
quantum pitch angle Oqu itself. For comparison the classical 
pitch angle is also shown. 




FIG. 7; The quantum pitch angle Oqu versus J2 for ferro- 
magnetic Ji and J± (case (iv)). The main panel shows the 
difference between the quantum pitch angle Qqu and its classi- 
cal counterpart Oci for various inter-chain couplings J± . The 
inset shows the corresponding data for the quantum pitch an- 
gle cvqu itself. For comparison the classical pitch angle is also 
shown. 



coupling {J± < 0.7 for case (i) and J± < 0.8 for case 
(ii)) we find a discontinuous behavior of the quantum 
pitch angle. A similar jumpwise change of the pitch an- 
gle has been found for 2D frustrated quantum spin mod- 
els [13, mi. Note that in Ref. [H the CCM curve was 
terminated before the jump occurs so that the jump has 
not been observed there [50|. The discontinuous change 
of aqu is related to the existence of two minima in the 
E versus a curve and takes place for parameter values 
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FIG. 8: The transition point J^ as a function of the inter- 
chain coupling J± for Ji = — 1 and J± > (case (iii)). 




FIG. 9: The quantum pitch angle Qqu versus J± for ferromag- 
netic Ji = — 1 (case (iii)) and various values of the frustrating 
NNN exchange J2. 



where both minima have equal depth, cf. Fig. [51 For 
larger J± (weaker quantum fluctuations) E{a) exhibits 
only one minimum and the quantum pitch angle is chang- 
ing continuously from the coUinear to the spiral GS. We 
remind the reader that for the classical case the transi- 
tion from the commensurate to the incommensurate GS 
takes place at J2 = |o^i|/4 and that it is continuous for 
any value of J±. 

Finally, in Fig. [5] we have drawn the pitch angle ofqu 
in dependence on the inter-chain coupling J± for case 
(i). As discussed above, with increasing J± the quan- 
tum pitch angle moves towards the corresponding classi- 
cal value. However, in accordance with the above discus- 
sion of the change of the sign of ofqu — ad we find two 
different regimes: For J2 < J2 the aqu increases with 
growing J± while for J2 > Jj the aqu decreases with J± . 

Let us now pass to cases (iii) and (iv), i.e. when 
Ji = —1 is ferromagnetic. By contrast to cases (i) and 
(ii), we find that here the transition from the coUinear 
state to the spiral state is always continuous. The rea- 
son for that can be again attributed to the strength of 



quantum fluctuations. For ferromagnetic Ji the rele- 
vant coUinear state at J2 < J2 consists of ferromagnetic 
chains, having a classical (i.e. 'non-fluctuating') GS, cou- 
pled by Jj^. While for case (iv) the ferromagnetic inter- 
chain coupling does not change this 'non-fluctuating' GS 
at all, for case (iii) due to the presence of antiferro- 
magnctic couplings quantum fluctuations become rele- 
vant and the GS becomes a true quantum state, but 
the change of magnetic correlations within the chains for 
Jj_ > remains weak. Hence, virtually no (case (iv)) or 
only weak (case (iii)) quantum fluctuations occur at the 
transition from the coUinear to the noncoUinear GS. As a 
result the transition takes place precisely at J2 = |<^i|/4 
for case (iv) while for case (iii) the transition point J2 is 
above the classical value indicating again an order from 
disorder effect. However, the shift of J2 is small because 
the quantum fluctuations are weak, see Figs. [6] and [S] 
Nevertheless, such a small shift of J2 due to a finite inter- 
chain coupling might be important for systems such as 
Li2ZrCu04 and Li2Cu02 [Oj llMl being near the transi- 
tion point. Note that Bader and Schilling [5l| first found 
that the transition point is fixed at J2 = |./i|/4, if the 
coUinear state is the classical ferromagnetic one. Note 
further that a similar change from a discontinuous to a 
continuous transition was discussed in Ref. [3^. 

The quantum pitch angle for cases (iii) and (iv) is 
shown in Figs. [H] and [T] Obviously, there is also a signi- 
ficant difference between the quantum and the classical 
pitch angle, however, it is smaller than for cases (i) and 
(ii). For Jj^ = the largest difference of about O.OQtt is 
found at J2 ~ 0.5. For case (iv) for all values of J^ < 
and J2 > 0.25 the qantum pitch angle is larger than the 
classical one. On the other hand, for case (iii) the shift 
of the transition point J2 leads to a change of the sign of 
«qu — ckci, i.e. the antiferromagnetic inter-chain coupling 
yields are more subtil change of the GS correlations by 
quantum fluctuations. Moreover, due to the shift of J 2 
we find a quite large difference jctqu — ctcil for J^ > and 
J2« 0.27... 0.28. 

Similarly as for cases (i) and (ii) we observe that for 
large J2 > Ji the model with strongest quantum fluctu- 
ations, i.e. Jj^ = 0, approaches the limit aqu — > 7r/2 
most rapidly. We mention here that a large value of 




for which 
The varia- 



1^2/^1! is realized e.g. in LiCuV04 

J2/J1 « —2.4 has been estimated 

tion of ckqu with J^ for various J2 is shown in Fig. [S] It 

can be seen that the variation of aqu with J±^ is largest 

for small J^. Furthermore, in difference to cases (i) and 

(ii) aqu decreases monotonously with Jj^ for all values of 

J 2 > Ji- 

As mentioned in Scct.[Tl a classical relation like Eq. ([T|) 
between the exchange couplings and the pitch angle has 
been used to discuss the J2/J1 ratio in Refs. [20, [2l| this 
way ignoring quantum effects. To overcome this problem 
we will provide an empirical formula that fits the contin- 
uous part of the aqu( J2) given in Figs. [U [21 [HIT] very well. 
Having in mind that the shape of the continuous part of 
the aqu(</2) resembles the classical behavior we find that 



TABLE I: Numerical values for the exponent v and the tran- 
sition point J2 in dependence on J± for the cases (i)-(iv), cf. 
Eq. HSl). 
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FIG. 10: The exponent v, see Eq. (|10p . in dependence on Jx 
for the cases (i)-(iv). 



ot<^\J\J'i) is well approximated by (written now in dimen- 
sional exchange units for the convenience of application 
in real experimental situations): 



*qu 



(J2) 



arccos 



'Ji 



[4( J2 - J| 



^il 



(10) 



with the exponent v as fitting parameter. Obviously, 
Eqs. (dni) and ^ coincide for J^ = | Ji|/4 and v = \. In 
Fig. [To] we show v in dependence on Jx for the four cases 
(i)-(iv). We find that v is always larger than the classical 
value Vc\ = 1. In accordance with the above discussion, v 
decreases with increasing Jj_, i.e. it goes towards the clas- 
sical exponent Vd- From the experimental point of view 
the edge-shared chain cuprates are of particular interest. 
The parameter situation of these compounds corresponds 
to case (iii). Hence, we give here simple fit formulas for 



that case which describe the behavior of J2{J±) as shown 
in Fig. [5] and the behavior of i'{J±) as shown in Fig. [5] 
by the red hne. We find that 



c^- 



yn 



+ bY 



J^ 



^..l..l 



tanh 



4i'<"' 



with a ^ 1.37, b = 0.13, c = 0.17, p = 0.036, and q ^ 2.11 
provides a reasonable fit of our data for case (iii). In ad- 
dition, we present numbers for J2 and v for various values 
of J± for all cases in table HI The fitting formulas pT|) 
as well as the data in table [I] can be used in combination 
with Eq. (fTU]) to fix the J2/J1 ratio using the pitch angle 
a as an input, e.g. from neutron scattering [9l. [lOl. [l3. |20| . 
Finally, we note that recently the pitch angle of 
Li2ZrCu04 has been determined from ^Li-NMR-data to 
amount a ~ 33° ± 2° [1^ . This value corresponds to a 
predicted ratio —J2/J1 = 0.298 within the framework of 
the classical Eq. (1). This value is surprisingly very close 
to the ratio —J2/Ji= 0.3 estimated from thermodynamic 
properties within the ID-quantum spin-1/2 Ji-J2-model 
[17| . However, from Fig. [5] for a realistic weak effective 
antifcrromagnctic inter-chain coupling J± < 0.1|Ji| a 
somewhat larger pitch angle would be predicted, namely 
a « 42° for —J2/Ji~ 0.3. Hence, other factors such as 
a sizable exchange anisotropy are expected to be rele- 
vant in this material. Such an anisotropy may lead to a 
modification of the classical spiral as well as to a reduc- 
tion of quantum fluctuation compared with the isotropic 
spin model considered here. The corresponding effects 
are outside the scope of the present paper. They will be 
considered elsewhere. 



IV. SUMMARY 

Based on the coupled cluster method (CCM) we have 
studied the GS correlations of a 2D array of frustrated 
spin- 2 J1-J2 chains coupled by an inter-chain exchange 
interaction Jx . Wc have discussed the influence of quan- 
tum fluctuations, frustration and inter-chain coupling on 
the pitch angle and the transition point between a GS 
with coUinear commensurate correlations and a GS with 
incommensurate spiral correlations. Using the CCM 
within the so-called SUB2-3 approximation we obtain a 
closed set of analytical equations which can be used to 
calculate the pitch angle for an arbitrary set of exchange 
parameters. We have found that for J2 > 0.25|Ji| the 
pitch angle of the quantum model significantly deviates 



from the classical value and can be strongly influenced 
by the inter-chain coupling J±. Furthermore, we have 
observed that the quantum pitch angle approaches its 
limiting value tt/2 for increasing values of J2 much faster 
than for the classical model. For several combinations 
of the sign of Ji and J± we have found that the change 
of the pitch angle within the quantum model could be 
discontinuous while the change of the pitch angle within 
the classical model is always continuous. 
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APPENDIX A: THE CCM SUB2-3 
APPROXIMATION 
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FIG. 11: Illustration of the five configurations on the spin 
lattice which contribute to the CCM SUB2-3 approxima- 
tion. Each configuration is related to a particular multi- 
configurational creation operators {C^} and to the corre- 
sponding correlation coefficient ol, L = 1,...,5, see Eqs. 
([5]) and ((8]). The circles in the figure represent lattice sites, 
the black circles indicate the position of the fiipped spins in 
a certain configuration. 



Within the SUB2-3 approximation scheme the CCM 
correlation operator S, see Eq. ([5]), contains five non- 
cquivalcnt correlation coefficients a^, L = 1, ... ,5, cor- 
responding to the lattice configurations shown in Fig. 1111 
These configurations (or lattice animals) represent the ar- 
rangement of spin operators acting on the lattice spins. 
To each configuration belongs a corresponding kct equa- 
tion, see Eq. ^. The set of these ket equations can be 
found by a bit tedious but straightforward calculation. 
It reads 



^{[cos(a) - 1](1 - 12a? + 8al + 8al + 8al + 24a'l) + 4a2[cos(a) + 1] - 8ai cos(a)} + -^{[cos(2a) - I](16a3a5 



Ji 

4 



jF jA 

-16aia2) +4ai[cos(2a) + 1] - 16aiCos(2a)} + ^(Sas - 8ai) ^(160405 - 16aia3 + 16a2a5 - 8ai) = (Al) 



— {[cos(a) — 1]( — 16aia2 + 160305) + 4ai[cos(a) + 1] — 16a2 cos(a)} H {[cos(2q) — 1](1 — 12a2 + 8ai 

jF jA 

+803 + 804 + 1605) — 8a2 cos(2a)} ^5^(802) i (160105 — 160203 — 802) = 



(A2) 



— {[cos(Qf) — 1]( — 160103 + I60205 + 160405) +8o5[cos(o) + 1] — 16o3Cos(q)} + — {[cos(2a) — 1]( — I602O3 

jF jA 

+ I60105) - I603 cos(2a)} + ^(4o4 - 403) - -^(1 - I203 + 80? + 802 + 804 + 24o5 - 403) = 



(A3) 



— {[cos(Qf) — 1]( — I601O4 + I603O5) — I604 cos(a)} H {[cos(2a) — 1]( — I602O4 + 805) 

jF jA 

— 16o4cos(2a)} H — ;^(4o3 — 804) ;^(— 160304 + I60105 — 804) = 



(A4) 



— {[cos(a) — l](16a203 — I601O5 + I603O4) + 803(008(0) + 1] — 32o5 cos(q)} H {[cos(2a) — l](16oi03 — 3202O5 

jF jA 

+ I60405) + 8o5[cos(2a) + 1] — 3205 cos(2a)} H ^(801 — I605) (160102 — I60305 + I601O4 — I605) — . (A5) 



r 



The GS energy is a function of (some of) these corre- lation coefficients and of the pitch angle a. It reads 
I 



E = ^{2ai[cos(a) - 1] + cos(a)} + ^{2a2[cos(2a) - 1] 



,(2a)]+'^--^{4a, + l) 



(A6) 



r 



Note that in the above equations one has to set J^ = 
(Jj^ = 0) and to replace J^ by J±_ [J^ by Jj_) for fer- 
romagnetic (antiferromagnetic) J^. To determine the 
quantum pitch angle ttqu as a function of the param- 
eters Ji, J2 and Jj^ one has to solve the equation 
dE / da\a=a u = 0, cf. Sect. [Ill together with the set 
of ket equations (JA1|1 - (jA5|l selfconsistently by standard 
numerics. 



Finally, we will illustrate some limiting cases contained 
in Eqs. (jAip - (|A6p . For Jj_ = one has 03 = 04 = 05 = 
and the remaining two non-trivial Eqs. (jAip and (|A2p 
then coincide with the corresponding equations given in 
Ref. [23. In case of J2 = 0, Jj^ = 0, 03 = 04 = 05 = 
but Ji 7^ (or alternatively, J2 = 0, Ji = 0, ai = 02 = 
as = but Jj^ 7^ 0) one finds the two ket equations for 
the simple unfrustrated linear chain. 
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